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ABSTRACT 

We derive the disk /-band luminosity function from the Zheng et al. sample of ~ 1400 disk M dwarfs observed 
with the Hubble Space Telescope. We adopt a Galactic-height-dependent color-magnitude relation to account 
for the metallicity gradient above the Galactic plane. The resultant /-band luminosity function peaks at Mj ^9.5 
and drops sharply toward Mj ~ 10.5. 

Subject headings: stars: late-type - stars: low-mass - stars: luminosity function - stars: statistics - surveys 



1. INTRODUCTION 

The stellar luminosity function (LF) measures the number 
density of stars as a function of luminosity. It is traditionally 
expressed in the V-band. However, sometimes it is more con- 
venient to use the /-band LF for planning and modeling ob- 
servations. Here we present the /-band LF of disk M dwarfs 
fr om Hubble Space T elescope (HST) star counts. 

IZheng et all J2001I) studied a sample of about 1400 disk M 
dwarfs found in 148 fields observed with the Wide Field Cam- 
era 2 (WFC2) on the HST and 162 fields observed with pre- 
repair Planetary Camera 1 (PCI). The V-band LF and the 
Galactic disk parameters were derived simultaneously using 
the method of maximum likelihood. If the luminosity of M 
dwarfs depended only on their color, as is the case for example 
for the solar neighborhood color-magnitude relation (CMR) 
determined bv iReidl lll991l) . the conversion from the V-band 
LF to the /-band LF would be straightforward. However, since 
many stars i n their sample lie far from the Galactic plane, 
Zhen g et all J2001I) introduced a more realistic CMR [CMR 
(2) in that paper], which makes a first-order correction for the 
metallicity effect: the luminosity then depends not only on the 
color but also on the height, z, above the Galactic plane. With 
this more complex CMR, the transformation from the V-band 
LF to the /-band LF is not trivial. 

In this paper, we apply the method of maximum likelihood 
to the same set of HST observations as in Zheng et al. (2001) 
and incorporate the Galactic-height-dependent CMR to derive 
the /-band LF of M dwarfs. This direct method avoids the un- 
certainty and the difficulty in the transformation from the V- 
band LF to the /-band LF. For details on the observations, the 
CMR, and the maximum likelihood method, we refer read- 
ers to lZheng et al.l J2001I) and references therein. In § 2, we 
present our results after a brief description of models. 

2. MODELS AND RESULTS 

The LF of disk M dwarfs is a function of the Galactic posi- 
tion (R, z). We parameterize it as the product of a vector of 12 
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absolute /-magnitude bins and a Galactic density law, 

R-R Q S 
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where $, is the solar neighborhood LF for the z'th magnitude 
bin, Ro = 8 kpc is the Galactocentric distance of the Sun, and 
H is the scale length of the disk. The vertical density profile 
v(z) is assumed to have either a "sech 2 " form, 



v,(z) = (1 - /?)sech 2 ±- + (3 exp ( -H 



or a "double-exponential" form, 



,<,!. ) = ( !- ;,exp(-Hj +/ 3exp(-g 
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The magnitude bins are centered at M/ = 6.5, 7.0, 7.5, 8.0, 
8.5, 9.0, 9.5, 10.0, 10.5, 11.0, 12.0, and 13.5, respectively. 
The size of each bin is 0.5 mag except the last two (1.5 mag). 
Given a choice of disk scale length H, the local LF, <P;, and 
the three disk profile parameters (h\, hi, (J) are found by maxi- 
mizing the likelihood. The scale length is determined by find- 
ing the maximum in the maximized likelihoods of an ensem- 
ble of solutions using different values of H. 

As in Zheng et al. (2001), we select stars below z= 2400 pc. 
The absolute /-magnitude of stars should satisfy 6.25 <Mj< 
14.25. The faint boundary is set to prevent contamination by 
spheroid giants. When translated int o color, the bound ary we 
use here is slightly bluer than that in Zheng et al. (200l|), but 
it is still acceptable since our results show virtually no change 
if we move the boundary to a higher absolute /-magnitude. 
Altogether, 1403 stars satisfy our selection criteria and they 
constitute our sampl e. This sam ple has 30 more stars than the 
one used by Zheng et al. ( 2001 ), which results primarily from 
the small difference in the blue boundary. 

We solve for the LF and disk profile parameters by max- 
imizing the likelihood and find that the resultant disk pro- 
file parameters, for either the sech 2 model or the double- 
exponential model , are almost identical to those derived in 
Zheng "eTall fcOOll) . This consistency of the results implies 
that the disk profile derived from the maximum likelihood 
method is robust. We therefore decide to fix the disk param- 
eters at the values in Zheng et al. (2001), namely hi = 270 
(156) pc, h 2 = 440 (439) pc, (3 = 56.5% (38.1%), and H = 2.75 
kpc for the sech 2 (double-exponential) model. 

Similar to the the case in V-band, due to the lack of local 
stars in the sample, the best-fit LF for the sech 2 model and that 
for the double-exponential model have nearly the same shape 
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Table 1. The /-band Luminosity Function 



M, <j> 
(mag) (ltrV^Z-mag" 1 ) 



6.5 


4.5±1.1 


7.0 


5.3±1.2 


7.5 


6.0±1.2 


8.0 


10.1±1.5 


8.5 


13.0±1.8 


9.0 


19.5±2.1 


9.5 


23.1±2.2 


10.0 


14.8±1.8 


10.5 


7.6±1.5 


11.0 


8.4±1.6 


12.0 


3.4±0.6 


13.5 


3.3±0.9 



but differ from each other in normalization, and the normal- 
ization difference is compensated by the vertical density pro- 
file. We adopt a linear combination of the LFs of these two 
models as our fin al result, with the s ame combination coef- 
ficient derived by Zhena et al. (2001) by matching the local 
normalization to the well-established luminous end of the so- 
lar neighborhood M-dwarf V-band LF. 

The /-band LF from our analysis is listed in Table 1 and 
shown in Figure 1. As expected, the overall shape of the I- 
band LF, with a peak at M/ ~ 9.5 and a sharp drop t oward 
Mj ~ 10.5, mimics the V-band LF in Zhen g et alJfcOOll) . 
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FIG. 1. — The /-band luminosity function based on 1403 M dwarfs observed 
with HST. 
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